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ABSTRACT 

We study the survival of ultrahigh energy nuclei injected in clusters of galaxies, as well as their 
secondary neutrino and photon emissions, using a complete numerical propagation method and a 
realistic modeling of the magnetic, baryonic and photonic backgrounds. It is found that the survival 
of heavy nuclei highly depends on the injection position and on the profile of the magnetic field. 
Taking into account the limited lifetime of the central source could also lead in some cases to the 
detection of a cosmic ray afterglow, temporally decorrelated from neutrino and gamma ray emissions. 

We calculate that the diffusive neutrino flux around 1 PeV coming from clusters of galaxies may 
have a chance to be detected by current instruments. The observation of single sources in neutrinos 
and in gamma rays produced by ultrahigh energy cosmic rays will be more difficult. Signals coming 
from lower energy cosmic rays {E < 1 PeV), if they exist, might however be detected by Fermi, for 
reasonable sets of parameters. 

Subject headings: ultrahigh energy cosmic rays: propagation, secondary emissions, clusters of galaxies 



1. INTRODUCTION 

The origin of ultrahigh energy cosmic rays still remains 
an unsolved puzzle. Most puzzling is the fact that at the 
highest energies, the cosmic ray arrival directions as ob- 
served by the current and previous experiments do not 
point back to any of their powerful candidate sources. 
One may underline that the Pierre Auger Observatory 
has reported a significant correlation between the arrival 
direction of cosmic rays of energy E > 5.7 x 10 19 eV and 
a catalogue of Acti ve Galactic Nuclei (A GN) located at 
less than 75 Mpc (| Abraham et all 12007b . This correla- 
tion however is not confirmed and should be interpreted 
as an eviden ce of an anisotropic d istribution of sources 
in the sky ([Abraham et al.l 120091 and see for exam- 
ple iKotera fe Lemoind I2008U iKashti fc Waxmanl 12 008. 
i Ghisellini et al.l J2008L ICuoco et al.l 120081 iTakami et all 
120081 iZawe t al. 2009Jfor interpretations of these results). 

Magnetic fields certainly play a key role in this puzzle. 
Yet our lack of knowledge concerning the distribution 
and strength of these fields, at both Galactic and extra- 
galactic scales, introduces a large degree of uncertainty 
in the subject of charged particle propagation. 

In view of this picture, one promising way to identify 
the sources of ultrahigh energy cosmic rays is to resort 
to multimessenger Astronomy. If charged particles are 
subject to deflections from magnetic fields, the direction 
of secondary photons and neutrinos 6 they produce by 
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interaction with the ambient medium should be much 
less affected during their propagation. 

Observing such secondary particles by using current 
detectors is not a straightforward task. Even though the 
existence of cosmogenic neutrinos is guaranteed, their 
dctcctability by present and future experiments is not. 
Indeed, their flux crucially depends on astrophysical 
parameters such as the cosmological evolution of the 
cosmic ray luminosity, the maxim um energy at the 
sources or the chem ical composition (j Allard et al.l 120061 : 
ITakami et al.l l2009). As for secondary high energy pho- 
tons, many studies have addresse d the question of their 
detec t ion dFerrigno et all 120041: iGabici fc Aharonianl 
120051 I2007t lArmengaud et al.l l2006h . but the situation 
remains unclear. These works indeed conclude that 
the observation of these photons strongly depends on 
assumptions on the configuration and strength of the 
intergalactic magnetic fields, as well as on the source 
luminosity in cosmic rays. 

Clusters of galaxies can be considered in various means 
as an ideal secondary particle production region. First 
they are dense regions of the Universe that harbour 
many candidate sources for ultrahigh energy cosmic rays, 
such as AGN, compact stellar remnants and cosmological 
shocks. We will discuss the importance of these possi- 
ble accelerators in our framework in section l2~2l Second, 
clusters of galaxies have enhanced photonic and baryonic 
backgrounds as compared to the extra-cluster medium. 
Indeed, X-ray observations of the bremsstrahlung emis- 
sion in clusters of galaxies have revealed that a large frac- 
tion of the cosmologica l baryonic density was contained 
in the intracluster gas (| White fc Fabianlll995h . Besides, 
the high concentration of ga laxies produces a loc al over- 
density of infrared photons (Lagache et al. 20*051). 

The most outstanding feature of clusters of galaxies 
for high energy cosmic ray propagation is certainly their 

duced by interactions of cosmic rays with the ambient medium as 
secondary photons and neutrinos. 
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strong magnetisation. Clusters of galaxies are the only 
extragalactic regions where magnetic fields could actu- 
ally be observed, by measuring the Faraday rotation of 
linearly polarised emission. They can reach values rang- 
ing from a few micr ogauss on scales of order ~ 10 kpc 
for normal c lusters (|Kim et al.lll991fc iClarke et alJl2001t 
IClarkdl200l . up t o 10 - 40 ulG on scales of 3 - 5 kpc for 
cool core clusters ([Taylor fc Perlevl[l993t lEnfilin fc Vo"gH 
2006). Such strong magnetic fields can easily confine 
cosmic rays of energy E < 5 x 10 17 eV x Z , where Z 
is the electric charge, for several hundreds of millions of 
years, leading galaxy clusters to act like storage rooms 
for these particles. While propagating inside the cluster, 
cosmic rays can thus experience many interactions with 
the enhanced photonic and baryonic backgrounds, and 
produce secondary neutrinos and photons that might be 
detectable with current and upcoming instruments. 

Many analytical and semi-analytical works were con- 
ducted on this subject to calculate the secondary gamma 
ray and neutrino fluxes, and to investigate the con- 
tribution of high energy protons produced in clusters 
of ga la xies to the total s pectrum dDar fe Shavivi 119951 
1996; iBerezinskv et al.l 119971: IColafrancesco fe Blasil 
19981 Ide Marco et all 120061: Urmengaud et alj 120061 : 
Murase et al.l I2008at I Wolfe et alJ 12008( 1 Estimates 



these fluxes can be calculated by assuming a proton in- 
jection spectrum, a Kolmogorov diffusion regime in the 
magnetic field, and a baryonic background profile. Nu- 
merical propagations of ultrahigh energy protons in more 
realistic three dimensional cl u ster m agnetic fields were 
also studied by iRordorf et all (|2004l ) . 

It can be summarised from these studies that it would 
be difficult to detect ^PeV energy neutrinos emitted 
from single sources with the upcoming generation of 
telescopes such as IceCube and KM3Net. The cumu- 
lative emission however may have some chance to be 
observed, depending on the source luminosity in cos- 
mic rays and on the background densities. GeV/TeV 
gamma ray signals may also be observed, again depend- 
ing o n assumptions on phy s ical p arameters: in partic- 
ular, IColafrancesco fc Blasil (|1998| ) find that the cluster 
population can account for a fraction ~ 0.5 — 2% of the 
diffuse gamma ray background and derive a list of clus- 
ters that could be obser vable in gamma r a y wit h the 
Fermi Space Telesco pe. lArmengaud et al.l ()2006l ) and 
iMurase et al.l (|2008af ) obtain less optimistic results, but 
the former authors find that a source located at 100 Mpc 
with maximum proton energy E max ~ 10 21 eV and in- 
jection spectral index ranging from 2.3 to 2.7 can be de- 
tected by HESS2. 

One should note that these previous studies generally 
assume a spherical symmetry for the modelling of the 
galax y cluster baryonic de nsity and magnetic field (apart 
from iRordorf et~atll20p who propagate particles in a 
simulated cluster magnetic field). The infrared photon 
background was taken into account only in the work by 
Ide Marco et~aTl ([2006), but in a highly optimistic way. 
Besides, all these studies investigate the propagation of 
protons inside clusters of galaxies. 

The chemical composition of high energy cosmic rays 
is still an open question though. The KASCADE data 
and measurements prior to the Pierre Auger Observatory 



indicate a dominance of heavy nuclei around the knee re- 
gion followed by a transitio n towards a light e r composi- 
tion around E ~ IP 18 5 eV (lApel et alj 12008 1 : iBird et all 
119931 : iShinozaki fc et aD 120051 : 1 Abbasi et al.ll2005[ ). The 
results of the Pierre Auger Observatory tend to suggest 
a mixed composition at all energies, that mig ht even get 
heavier at the highest end (jUnger et al.ll2007l ). 

The energy losses of nuclei with mass number A > 1 
during the propagation definitely differs fr om those of 
proton s. They were first studied in d etail bylPuget et al.l 
(1976) a nd have been reexamin ed by lStecker fc Salamonl 
Il999li. lEpele fc Rouleti (| 19981 ) and most recently by 
IKhan et al.l (|2005[ 7! The implementation of these re- 
sults analytically or in detailed numerical simulations 
for the propag ation of ultrahigh energ y cosm ic rays 
(jBertone et ali l2002T: lAllard et al.l I2005L 120(161 . I2007bf > 
have demonstrated that iron nuclei can travel hundreds 
of megaparsecs before losing their energy. This inter- 
esting property combined moreover with the fact that 
such nuclei can be accelerated to an energy typically 
Z times larger than protons, has motivated several au- 
thors to consider heavy nuclei as primary co s mic ra ys 
(lAnchordoaui et all 119991: lAllard et all 120061. I2007allbl: 
Tnoue et al.ll2007t lAnchordoaui et al.ll2008b lAllard et all 
20081) . lAllard et alj (|2008D conclude that whatever the 



composition injected at the source, the propagated spec- 
trum detected at the Earth should be either dominated 
by protons, or by heavy nuclei of iron type. The latter 
case occurs if the source is strongly enriched in iron or if 
the proton maximal injection energy is smaller than the 
GZK cut-off energy, so that only heavy nuclei are present 
at the greatest energies. 

It is thus of prime importance to have an idea of the 
composition that escapes from the sources of ultrahigh 
energy cosmic rays, and this question should also be 
addressed in the context of propagation inside clusters 
of galaxies, considering that they might host powerful 
accelerators. There is no reliable prediction of the 
expected composition at the source, mainly because 
very little is known on the physical parameters that 
govern the acceleration and survival of nuclei in those 
powerful objects. As a case study nonetheless, one 
might reasonably consider a composition similar to the 
Galactic cosmic ray one, that contain around 30% of 
CNO and heavier nuclei. Moreover, AGN that we will 
consider as sources of cosmic rays in our context, are 
known to have supersolar m etallicities (jGroves et alJ 
12001 iMathur fc Fields! l200l ). which further justifies 
our choice of injecting a mixed composition. Because 
the diffusion time in a magnetic field increases with the 
charge Ze of the particle, heavy nuclei should remain 
confined longer times in the structure, leading to an 
enhanced number of interactions and thus possibly to 
complete depletion of the original particle. Distinct 
signatures of these propagation effects could be observed 
in the produced spectrum as well as in the neutrino and 
gamma ray fluxes. 

In the present paper, we explore the consequences of 
the injection of a mixed chemical composition in galaxy 
clusters, by calculating the propagated primary and sec- 
ondary particle fluxes. We elaborated for this purpose a 
propagation code that treats all baryonic and photonic 
interaction processes for primary and secondary nuclei, 
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as well as their accurate trajectories in the magnetic field, 
gamma ray signals formed by electromagnetic cascades 
of ultrahigh energy pairs and photons in the intergalactic 
space are calculated in a second step as a post- analysis. 
We model clusters of galaxies using three dimensional 
output s of magneto hydrodynam ical (MHD) simulations 
run by iDubois fe Tevssierl (|2008l ) , making a distinction 
between cool core and non cool core clusters. 

This paper is organized as follows: in section [2l we 
will describe our physical modelling of galaxy clusters for 
high energy cosmic ray propagation. The details of the 
numerical techniques of our propagation code are given 
in section [3) We will present and discuss our results are 
presented in section [4j 

2. MODELING CLUSTERS OF GALAXIES 
FOR ULTRAHIGH ENERGY COSMIC RAY PROPAGATION 

We discuss in this section our modelling of clusters 
of galaxies for ultrahigh energy cosmic ray propagation. 
This comprises a three dimensional modelling of the mag- 
netic field and infrared photon background, a consistent 
baryonic background profile and an adequate choice of 
sources for injection. Features for both cool core and 
non cool core clusters will be discussed. 
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Fig. 1. — Magnet i c field intensity profiles from the simulations of 
IDubois & Teyssier (20080. The red solid line represents the mean 
magnetic field intensity in radial bins, for the cool core cluster of 
galaxies. The blue dashed line is for the non cool core cluster. The 
magnetic field intensity is rescaled for our simulations to obtain a 
field of 1, 3 or 10 fiG at the center for cool core clusters, and 1 fiG 
for non cool core clusters, see text for details. 



2.1. Magnetic field 

Clusters of galaxies can be roughly split into two cat- 
egories: cool core clusters and non cool core clusters. 
Observations indeed indicate that half of the cluster pop- 
ulation exhibits a peaked X-ray emission at the core and 
a central temperature of 30 to 40% of the y irial tem- 
perature (jChen et al.ll2007t iBurns et al.H "2008). This bi- 
modality has also some significant consequences when it 
comes to modelling magnetic fields: cool core clusters are 
reported to have a much stronger magnetic field in the 




core, as well as a higher tur bulence rate, as comp ared 
to non cool core clusters (see iCarilli fe Tavlori 2002 for a 
recent review). Such differences can play a non negligi- 
ble role in cosmic ray confinement, and necessitates some 
investigations. 

We model our cluster magnetic field using the three 
dimensional outpu t s of the MHD simulations run by 
IDubois fc Tevssierl (|2008l ). The simulations were run in- 
cluding Dark Matter, gas, ultraviolet heating, hydrogen 
and helium cooling, star formation and magnetic fields 
with the Adap tative Mesh Refinement code RAMSES 
(|Tevssierl [2002) . The galaxy cluster has been evolved 
up to z — within a 'standard' ACDM cosmology 
with parameters fl m — 0.3, £1a = 0.7, £1& = 0.04, 
Ho = 70kms _1 Mpc _1 and erg = 0.9. The cool core 
cluster has the following properties at z = 0: i?2oo 
(the radius inside which the mean interior overdensity 
is equal to 200) equals l.l/i Mpc, the virial mass 
M 20 o = 200 x (47r/3)p c i?| 00 = 3.5 x 10 14 /i^M© (where 
p c is the critical density), and the X-ray temperature 
T x = 5.1keV. 

Though this simulation enables one to have a very ac- 
curate resolution of 1.2 /i -1 kpc at the center of the clus- 
ter, the outputs were smoothed for practical reasons on 
a cubic grid of 256 3 regular cells and of size 5/i _1 Mpc, 
leading to a constant 19.5/i _1 kpc resolution. The au- 
thors performed two simulations: one including atomic 
cooling and star formation, which led to the formation of 
a cool core cluster, and the other without these processes, 
leading to a non cool core cluster. As a result of the in- 
crease of density due to the cooling flow, the magnetic 
field is stronger in the cooling run than in the adiabatic 
run by more than one order of magnitude, reaching a typ- 
ical value of 0.5 /iG in the most favourable case (see the 
magnetic field radial profiles in Fig. [1]). In both cases, 
compression and shear motions of the the hot plasma 
medium amplify and give rise to a non-trivial topology 
of the field lines at large scales (see Fig. [2]). Thus, it is of 
particular interest to propagate particles in these realis- 
tic magnetic fields that can strongly vary on kpc scales 
(see the typical coherence lengths on Tab. [I}. 

Many studies indicate that AGN jets can grow into 
large cavities through the hot plasma of the ICM 
(jArnaud et al.lll984f ). iBinnev fc Tabor] (| 19951) show that 
those high velocity jets could disrupt the formation of 
cooling flows in the core. This is a popular explana- 
tion to the non detection by Chandra and XMM of line 
emissions of cold gas that would have been expected in 
a standard cooling flow picture at cluster cores. Such 
re-heating process e s are not included in the work of 
IDubois fc Tevs sier (2008), but drawing a complete pic- 
ture of a self-consistently growing supermassive black- 
hole powering a high-velocity jet including magnetic 
fi elds still remains a c hallen ge in numerical treatment. 

IDubois fc Tevssierl (|2008f) scale t heir magnet i c fiel d 
seeds according to previous works bv lDolag et al.1 (2005): 
the initial comoving magnetic field Bq = B(z)/(1 + z) 2 
was set up to an arbitrary low value and then renormal- 
ized to 10 _11 G. It can be noted however that these au- 
thors consider non cool core clusters. It is thus justified 
to perform a new normalization of the field in adequacy 
to the type of the cluster that we study - up to the point 
where magnetic fields become dominant in the fluid dy- 
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Fig. 2. — From left to right: magnetic field intensity, infrared photon density and baryon density projected along the observed plane, 
from the simulation outputs of Dubois & Tcvssier (2008), in the case of a cool core cluster. The colour contours are in logarithmic scale 
and the box size is of ~ 7 Mpc. The red line on the magnetic field slice represents the trajectory of one proton of energy E = 10 17 eV 
propagated in our simulations. The particle reaches the edge of the box in a time ~ 200 Myr. 



namics, which is not the case with the values used here. 
For cool core clusters, we normalise the maximum in- 
tensity of t he field to 3, 10 and 30 /iG (following values 
derived bv lEnfilin felfogtl |200H and we take a value of 
1 /iG at the center of non cool core clusters, that corre- 
sponds to the same scaling as for the cool core case at 
30 /iG (see figure [1]). 

The coherence lengths of the magnetic field, which are 
of prime importance for cosmic ray propagation, were 
calculated directly from the simulated fields, though 
assuming spherical symmetry for simplicity. We split 
the cluster in five radial shells and assume that in each 
of these shells, the coherence length is roughly constant. 
The magnetic field at a given distance r from the 
source can be decomposed into a global and a turbulent 
components: B(r) = (B)(r) + SB(r). We calculate the 
magnetic turbulence energy spectrum in each radial 
shell by applying a Fourier transform to the ratio 
B/(B) = 1 + SB /(B). We assume in this calculation 
that the purely turbulent component is defined by an 
isotropic power law spectrum. This method allows 
us to get rid of the linear dependance of SB(r) over 
(B)(r) and draw a spectrum of the pure turbulence 
energy. The coherence length is then calculated using 
the following formula (which is valid as long as the 
turbulent spectrum is steep enough): A ~ 0.77/ A; m i n , 
where fc m ; n is the wave number at which the energy 
spectrum reaches its maximal valu e, setting in the 
turbulent cascade IjCasse et alj 12002). The propagation 
method implemented in our code enables one to take 
into account turbulent effects of scale inferior to the grid 
resolution of the simulation (more details in section [3]) . 
We make use of this advantage by adding by hand a 
value for the coherence length in a shell of radius smaller 
than the resolution of the grid (r = 20 kpc). We set 
the value of this coherence length to 5 and 15 kpc for 
our cool core and non cool core clusters respectively, 
which are ave r age y alues according to the studies of 
lEnfilin fe Vogtl (|2006f) . Table Q] presents our calculated 
and assumed coherence lengths in each radial shell, for 
cool core clusters Acc and for non cool core clusters 
Ancc- 



TABLE 1 

Magnetic field coherence lengths in kpc 



min 


'"max 




Ancc 





20 


5 


15 


20 


100 


36 


10 


100 


200 


7.3 


80 


200 


800 


109 


109 


800 


7000 


145 


145 



Note. — Coherence lengths in each radial shell (delimited by 
the distances from the cluster center r m i n and r max ), for cool core 
clusters Acc an d for non cool core clusters Ancc- Values are in 
kiloparsecs. 



2.2. Possible ultrahigh energy cosmic ray sources in 
galaxy clusters 

As mentioned previously, clusters of galaxies are dense 
regions of the Universe that may harbour various can- 
didate sources of ultrahigh energy cosmic rays. Among 
those candidates, we are most interested in the frame- 
work of this paper, in those located in the inner dense 
regions of the cluster, where injected particles are bet- 
ter confined and have higher interaction rates, thus some 
chance of leaving observable signatures. 

Long duration gamma ray bursts (jWaxmanl 119951: 
lVietrilll995l : iMurase et all 120061 l2008bfl and s ome other 
compact objects like magnetars (|Aronsl I2003D are good 
candidate sources for ultrahigh energy cosmic rays. 
These transient sources are expected to be associated 
with the death of young massive stars which are com- 
mon in star-formin g galaxies (including dwarf and ir- 
regul ar galaxies, see lLe Floc'h et aLll2006t iSavaglio et al.l 
2009). Galaxies with high star formation rate are how- 
ever not predominant in galaxy clusters, and are partic- 
ularly rare in the central de nse part which is of interest 
for cosmic ray confi ne ment (iKoopmann fe Kennevl 120041 : 
IGavazzi et all 120061 ). IGavazzi et al I (|2006l) re"port that 
the star formation rate per unit mass of high luminosity 
spiral galaxies that are projected within 1 virial radius is 
about a factor of two lower than at larger clustercentric 
projected distances. 

Accretion and merger shocks in galaxy clusters were 
considered by many authors as i mportant cosmic ray ac- 
celerators (iNorman et all 119951: [Kang- et al.l 119961 1 1997 : 
iMiniati et al.ll2000t iRvufcK ang 20031; llnoueet al.ll2005 . 
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120071 iMurase et al.ll2008af) . However, high Mach num- 
ber shocks that are needed to accelerate particles with 
a relevant spectrum mainly occur at the rim of clusters 
of galaxies, where the infrared photon and baryon back- 
grounds are very faint. Therefore, cosmic rays acceler- 
ated at such external shocks will experience less interac- 
tion s in clusters of ga laxies than if produced at the cen- 
ter (|Inoue et al.l 120071 ). In view of this situation, we will 
focus in this paper on AGN that are more commonly lo- 
cated in the interesting central region of clusters of galax- 
ies. 

What exact part and type of AGN are most suscepti- 
ble to accelerate particles, and which are more present in 
galaxy clusters are questions that are beyond the scope of 
this paper. We might underline however that radio obser- 
vations indicate a strong presence of radio-loud Faranoff 
Riley type I (FRI) galaxies in th e center of galaxy clusters 
(|Bestll200l : lGilmour et a l. 2009). In cool core clusters es- 
pecially, as mentioned earlier, these AGN are believed to 
shut off the cooling flow by mechanical energy dissipa- 
tion, which would explain the absence of cold gas emis- 
sion lines in the observations of Chandra and XMM. Pow- 
erful Faranoff Riley type II galaxies that are even better 
candidates for particle acceleration are more likely found 
in ga l axy groups and weak clusters ( Prestage fc Peacockl 
U988t Mill fc LillvHl991h iMiller et al.ll2002D . 

AGN of more quiet types (Seyfert galaxies for exam- 
ple) are also f ound at an average rate of 1.5 AGN per 
galax y cluster (jGilmour et al.l 120071 . 120091 : iMartini et all 
2006). These studies indicate that these common AGN 
are predominantly detected in areas of moderate density. 
For this reason, and mostly because such AGN are much 
less efficient to accelerate particles at very high energy, 
these quiet objects are of less interest in our framework. 
One might note nevertheless that they may still acceler- 
ate protons and nuclei up to high energies, for example 
in the vicinity of the central black hole or in weak jets. 
Moreover, it might be possible that some of such AGN 
were more active in the past, and produced cosmic rays 
that are still stored in the cluster. These particles could 
also act as see ds for further accelera tion at accretion and 
merger shocks ([Murase et alJl2008a| ). All these situations 
introduce many unknown parameters and are not in the 
scope of this paper. 

2.3. Photonic and baryonic backgrounds 

We consider in this study the interaction of high energy 
protons and nuclei with the Cosmic Microwave Back- 
ground (CMB) and with the infrared background - the 
term 'infrared' here includes as well ultraviolet and opti- 
cal backgrounds. We calculate the infrared background 
by adding the contribution of the diffuse extracluster in- 
frared photon density and the intracluster density cre- 
ated by galaxies inside the cluster. We model the dif- 
fuse infrared backgr ound according to the studies of 
iKneiske et ail (|2004f) and IStecker et alj (|2006fl . In this 
work, we do not include redshift evolution, as its effect 
is negligible as compared to the uncertainties on all our 
other parameters, especially as we consider only redshifts 
ofz<0.2. 

The intracluster infrared background is modelled as 
follows. We assume that local (z < 0.2) clusters of 
galaxies are mostly populated by elliptical galaxies that 
moderately enrich the cluster with infrared photons. We 



thus scale an average elliptical galaxy SED (Spectral En- 
ergy Distribution) with the galaxy density in our simu- 
lation cubes, convoluted by a 1/r 2 function in order to 
account for the flux diminution with distance. The SED 
was kindly provided by He rve Dole and details of their 
modelling can be found in lLagache et alj (|2005f ) and in 
references therein. The colour contours of the projected 
infrared photon density are drawn in the middle panel of 
Fig.H 

It is important to stress again that the galaxies 
in galaxy clusters are ellip ticals in maj o rity. Ac- 
cording to the a n alysis of iHolden et al.l (|2007l ) and 
Ivan der Wei et all (|2007l ). the fraction of ellipticals 
among galaxies of mass M > 4 x 1O 1O M in clusters 
at redshift z < 0.2 is higher than 90%. These galax- 
ies have SED that highly differ from those of ultra- 
luminous infrar e d gal axies (ULIRG) that were used by 
Ide Marco et al.l (|2006f ) to evaluate their cluster infrared 
background. 

Another important process for cosmic rays propagat- 
ing in the intracluster medium, especially at energies 
E < 10 18 eV, is their interaction with the cluster gas. 
The baryonic bac kground is evaluated usin g the simula- 
tion outputs from lDubois fc Tey ssicr lj2008f l. In order to 
benefit from the high resolution of the simulation in the 
densest region of the cluster, we do not use a smoothed 
out three dimensional mapping of the density grid, but 
the mean radial profile instead, with high resolution at 
the center. We checked that the density distribution is 
actually smooth enough over the whole simulation box, 
so that the loss of small features while averaging over 
radial bins does not affect our results. 
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Fig. 3. — Baryonic d ensity profiles from the simulations of 
IDubois &; Tcyssier (2005). The red solid line represents the mean 
baryonic density in radial bins for the cool core cluster of galaxies. 
The blue dashed line is for the non cool core cluster. 



3. PROPAGATION OF PROTONS AND NUCLEI: 
NUMERICAL TECHNIQUES 

We describe in this section some aspects of our ultra- 
high energy nuclei propagation code, that combines a 
fast and accurate semi-analytical trajectory integration 
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Fig. 4. — Mean free paths as a function of cosmic ray energy for various interactions with various backgrounds, all calculated in the 
inner 20 kpc region of the cool core cluster. Left panel: values for a propagating proton. Pion photo-production on the CMB (black solid 
line), on the diffuse extragalactic infrared background (red thin dotted line), on the cluster infrared background (red thick dotted line) 
and hadronic interactions on the baryonic gas (blue dashed line). Middle panel: comparison between pion photo-production process (black 
solid) and photo-disintegration for iron (red dashed) for interactions with the cluster infrared background. Right panel: variation of the 
hadronic interaction mean free paths for (from top to bottom) proton, helium, carbon, silicon and iron. The cross sections used in these 
calculations were provided by EPOS. 



method in the magnetic field, and complete Monte Carlo 
calculations of photonic and baryonic energy losses for 
primary and secondary nuclei. 

The transpo rt scheme in the magnet ic field was 
adapted from IKotera fc Lemoini (|2008aD . In this 
method, particles propagate by crossing spheres of 
diameter the local coherence length of the magnetic 
field. The exiting time and deflection angle of the 
particle when leaving the sphere are randomly computed 
following analytical distribution functions depending on 
the field strength, coherence length and Larmor radius 
of the particle. These functions include information on 
the small scale turbulence component via the diffusion 
co efficient that we im plemented according to the results 
of iCasse et al.l (|2002D . The advantages of this method 
are twofold: the trajectory integration is very fast as 
compared to a direct integration, and second, it allows 
to enlarge artificially the range of scales on which the 
magnetic field is distributed, and thus to account for 
turbulence scales below the simulation grid resolution. 
More det ails on this scheme can b e found in the ap- 
pendix of IKotera fe Lemoinel (|2008afl . 

The interactions of protons and nuclei with CMB, 
infrared, optical and ultraviolet photon backgrounds 
were mostly modelled ac c ordin g to t he Monte Carlo 
methods of lAllard et all (|2005l . I2006f) . These inter- 
actions produce features in the propagated ultrahigh 
energy cosmi c ray spectrum such as the "GZK cut- 
off" ljGreisenlll96a iZatsepin fe Kuzminl[l96l and their 
decay products generate the cosmogenic neutrino flux 
(jBerezinsky fe Z atsepin Hl969h . During the propagation 
in the intracluster medium, the same interactions are also 
expected to take place, at a rate depending on the am- 
bient photon density and the confinement time in the 
densest regions, that were described in the previous sec- 
tion. These interactions can also possibly produce spec- 
tral and composition features in the cosmic rays as well 
as secondary neutrinos and photons. 



In our calculations, electron and positron pair pro- 
duction from the interactions with the various photon 
backgrounds (also known as the Bethe-Heitler process) 
are treated as a continuous energy loss process. We use 
the cross sections, inelasticities, as well as the mass and 
charge scaling o f the at t enuat ion length for nuclei that 
are provided by IRachenl (|1996D . 

If the energy of the background photons exceeds 
~145 MeV in the nucleon rest frame, protons and neu- 
trons can interact through the pion photoproduction pro- 
cess. This latter process is taken into account in our code 
using Monte Carlo calcula tions based on the outputs of 
event generator SOPHIA (jMiicke et al.l 12000). The use 
of SOPHIA allows us to treat accurately the various in- 
teraction channels (direct pion production, resonances, 
multi-pion pr oduction) and their branching ratios (see 
IRachenl [19961) . 

The interactions experienced by nuclei with photon 
backgrounds differ from those of protons. In addition to 
the pair production losses that result in a decrease of the 
Lorentz factor and the rigidity of the nucleus, one must 
consider the photo-disintegration processes that lead to 
the ejection of one or several nucleons from the nucleus. 

Different photo-disintegration processes become dom- 
inant in the total interacti on cross section at different 
energies (|Puget et al.lll976| ). The lowest energy disinte- 
gration process is the Giant Dipole Resonance (GDR) 
which results in the emission of one or two nucleons and 
a particles. The GDR process is the most relevant as it 
has the highest cross section and the lowest thresholds, 
between 10 and 20 MeV for all nuclei. For nuclei 
with mass A > 9, we use the the oretically ca l culate d 
GDR cross sections presented by iKhan et al.l |2005), 
which take into account all the individual reaction 
channels and are in better agreement with data than 
previous treatments. For nuclei with A < 9, we use 
the phenomen ological fits to the data provided by 
IRachenl ([1996) ■ Around 30 MeV in the nucleus rest 
frame, and up to the photopion production threshold, 
the quasi-deuteron (QD) process becomes comparable 
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to the GDR and dominates the total cross section at 
higher energies. The photopion production (or baryonic 
resonances, BR) of nuclei becomes relevant above 150 
MeV in the nucleus rest frame (e.g., ~ 5 x 10 21 eV 
in the lab frame for iron nuclei interacting with the 
CMB). We use the parameterisation given by iRachenl 
(1996) in which the cross section in this energy range 
is proportional to the mass of the nucleus. For more 
details on the num erical method, one might refer to 
lAllard et alJ ((20051) . 

In addition to interactions with soft photon back- 
grounds, protons and nuclei can experience hadronic 
interactions with the baryonic background in the dense 
regions of the intracluster medium. We implement this 
process using the h adronic interaction model EPOS 
(jWerner et all [2009) . We assume that the baryonic 
background is essentially made of protons at rest, which 
is a good approximation in view of the extremely high 
Lorentz factors of the projectile nuclei. Our choice of 
using the new interaction model EPOS was motivated 
by the fact that it is the only model for very high energy 
cosmic ray interactions which is actually also used for 
particle physics analysis. It has indeed been tested 
thoroughly against most of the existing hadron-hadron 
and hadron-nucleus data. As a consequence, EPOS has 
a very good description of secondary products, which 
are of interest to our work. In EPOS, the fragmentation 
of heavy nuclei is not treated: spectator nucleons are 
set as free particles and the building up of resulting 
nuclei has to be done as post-processing. We make 
use for this purpose of the fragmentation model of 
ICampi &; Huefnen (|1981h . as i mplemented in the ai r 
shower simulation code CONEX (Bcrgm ann et al.ll2007f ). 

Figure |4] shows the mean free paths for the interactions 
described above, as a function of the cosmic ray energy, 
for the cool core cluster. Values for the non cool core 
case can be derived from these plots by simple scaling 
laws. 

The first plot presents the mean free paths for inter- 
actions experienced by protons, namely pion production 
and hadronic interactions - mean free paths for pair pro- 
duction are not represented. We draw separately the 
contribution of the CMB photon background (black solid 
line), the diffuse extragalactic infrared background (thin 
dotted red line) and in thick dotted lines the infrared 
background of the cluster itself, in the inner region (with 
distance from the center r < 20 kpc). It clearly appears 
that the interactions with the CMB photons dominate 
at extreme energies (E > 10 20 eV). One can also notice 
that the infrared background at the center of the cluster 
is highly enhanced as compared to the diffuse extragalac- 
tic background, resulting in a shortening of the mean free 
paths around E ~ 10 18 eV of more than 4 orders of mag- 
nitudes. One should however note that the extragalactic 
background becomes dominant for energies E < 10 19 eV 
at ~ 1 Mpc of distance from the center of the cluster. 
Obviously, the contribution of hadronic interactions will 
be predominant in the center of the cluster, at all en- 
ergies. Mean free paths for proton-proton interactions 
at different distances from the center can be easily ex- 
trapolated from this plot by scaling with the baryonic 



density presented in Fig. [3J Due to the rapid decrease 
of the density with the distance from the center of the 
cluster, hadronic interactions become quickly negligible 
and are completely overcome by photonic interactions at 
0.1 Mpc. 

The figure in the middle panel compares the mean free 
paths for interactions of protons (black solid line) and 
iron (red dashed line) with the cluster infrared back- 
ground, as a function of their energy. The distances 
are calculated again in the inner part of the galaxy 
cluster (r < 20 kpc). The mean free path for photo- 
disintegration of iron - that cumulates all GDR, QD and 
BR processes - are shorter by more than two orders of 
magnitude than for the pion photo-production with pro- 
tons. This implies that heavy nuclei can be strongly de- 
pleted by this process at the center of the cluster, but 
not necessarily that secondary photons and neutrinos 
will be emitted. Indeed, the only process here that pro- 
duces pions is the baryonic resonance (BR) which takes 
place only at extremely high energy. In comparison, pion 
production from proton-photon interactions can happen 
at all energies above E ~ 10 17 eV. The depleted nuclei 
can still produce secondary protons that might interact 
through this channel. It appears that accurate calcula- 
tions are definitely necessary to investigate the possible 
signatures of these different processes. 

The last panel of Fig. [4] presents the variation of the 
hadronic interaction mean free paths for a set of nuclei, 
as a function of their energy, in the inner 20 kpc of the 
cluster. The cross sections used in these calculations 
were provided by EPOS. It appears that at high energy, 
the mean free paths for nuclei can be derived from the 
proton-proton interaction case by a simple scaling law, 
but there is a slight discrepancy at lower energy. As 
expected, heavy nuclei have higher probability (of order 
A 2 / 3 times more, where A is the nuclei mass number) to 
interact with the cluster gas, leading to the production 
of a wide bunch of secondary particles. 

4. RESULTS AND DISCUSSION 

In this section, we present the results of our simula- 
tions for various normalizations of the magnetic field, in 
the case of cool core and non cool core clusters of galax- 
ies. We position the source of ultrahigh energy cosmic 
rays at the center of the cluster, or at a slightly shifted 
position, and inject either a 100% proton comp osition, 
or a mixed composition as in lAllard et all (|2006f L based 
on Galactic cosmic ray abundances. The injection spec- 
tral index is chosen to be 2.3 for all our plots (if not 
indicated otherwise) and a maximum injection energy 
of -E max ,z = Z x 10 20 ' 5 eV is assumed for a nucleus of 
charge number Z (an exponential cut-off is taken above). 
Most of our results are arbitrarily normalised to unity at 
E = 10 19 eV. When a comparison with observational 
data or limits are needed, we normalise our fluxes by 
choosing the cosmic ray luminosity of the central AGN 
to be L CT — 10 45 erg s -1 with a minimum injection en- 
ergy E m i n = 1 GeV. This integrated luminosity indeed 
allows to match the ultrahigh energy cosmic ray flux as- 
suming a source density of n a — 10~ 5 Mpc -3 , for this 
value of the spectral index. Moreover such a luminosity 
is reasonable for this type of object. 

4.1. Cosmic ray spectra 
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Fig. 5. — Cosmic ray energy spectra for central injection and 
for our cool core cluster with magnetic field at the center: B c = 
3 fiG, 10 fiG and 30 /iG from top to bottom. The black solid line 
presents the total cosmic ray flux and the blue dashed line the flux 
of secondary nuclei produced during the propagation. We also show 
the contribution of the different chemical species to the total flux 
as indicated in the second panel: protons (green), Helium (black 
dot-dashed), Carbon, Nitrogen and Oxygen (pink), all species with 
charge Z comprised between 12 and 20 (red) and heavy nuclei with 
charge 20 < Z < 26 (brown). These spectra are normalised to 



unity at E = 10 la eV. 



One key point of this paper is to investigate the result- 
ing composition of ultrahigh energy cosmic rays, when a 
mixed chemical composition is injected at the center - or 
at a slightly shifted position - of a cluster of galaxies. 

Figure [5] shows our resulting cosmic ray spectra for our 
cool core cluster of galaxies for various normalizations of 
the magnetic fields (at the center: B c = 3 fxG, 10 /J.G and 
30 fiG from top to bottom). The black solid line presents 
the total cosmic ray flux and the blue dashed line the 
flux of secondary nuclei produced during the propaga- 
tion. We also present the contribution of different chem- 
ical species to the total flux as indicated in the second 
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FiG. 6. — Same as Fig. [5] but for our non cool core cluster with 
B c = 1 fiG and central injection (top), and for our cool core cluster 
with B c = 10 fiG and the source positioned at 100 kpc from the 
center (bottom). 



panel: protons (green), Helium (black dot-dashed), Car- 
bon, Nitrogen and Oxygen (pink) , all species with charge 
Z comprised between 12 and 20 (red) and heavy nuclei 
with charge 20 < Z < 26 (brown). These three panels 
assume that the injection of particles occurs at the very 
center of the cluster, where the background densities are 
the highest. The limited lifetime of AGN is not con- 
sidered in these plots: we assume that the source emits 
particles continuously and that the stationary regime has 
been reached. 

In Figure [5] one can observe features in the spectra of 
the different nuclei, that increase in amplitude with the 
mass of the nuclei and the normalization of the magnetic 
field. At the lowest energies the flux of the heaviest nu- 
clei is greatly suppressed due to hadronic interactions. In 
the same energy range, light and intermediate nuclei suf- 
fer less interactions due to their larger mean free paths 
(Ia p oc A~ 2 / 3 ) and rigidities (as the confinement time 
Tconf oc (-E/Z) -1 / 3 ). As the baryonic density rapidly de- 
creases with the distance to the center and the confine- 
ment time with the energy, hadronic interactions become 
less efficient to disintegrate nuclei above E ~ 10 17 eV. 
Above this energy however, heavy nuclei start to inter- 
act with optical and near- infrared photons. The photon 
density decreases more slowly than the baryonic density; 
photo-disintegration is then very efficient to deplete the 
heavy nuclei component especially for the largest mag- 
netic field normalizations. For lighter nuclei, the energy 
threshold for photo-interactions is lower (more or less 
proportional to the mass of the nuclei A) but the mean 
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Fig. 7. — Evolution of the cosmic ray spectrum in time, assuming a lifetime of £agn = 10 Myr for the central AGN, for the case of a cool 
core cluster of central magnetic field B c = 10 fiG. Each panel presents the spectrum at the time indicated at the top-right hand corner. 
The injection from the source (AGN) is assumed to begin at t = 0. The contribution of the different chemical components are shown in 
the same color code as for Fig. \E\ The thick black line is the total spectrum and the thin black line indicates the total flux obtained for an 
infinite AGN lifetime and an integration of the flux over a Hubble time (stationary regime), as in Fig. [5] The spectra are normalised to 
the value of the stationary flux obtained at E = 10 19 eV. 



free path is larger (more or less proportional to A^ 1 ), 
therefore the effect of photodisintegration on light and 
intermediate nuclei is weaker. 

The production of secondary protons and their relative 
abundance at low energies depend on the interactions of 
nuclei. One can see as expected, that the relative abun- 
dance of secondary protons is higher for the largest mag- 
netic field normalization due to the better confinement 
in the vicinity of the cluster center. Between E ~ 10 175 
and 10 18 ' 5 eV the decrease of the secondary proton abun- 
dance results from the deconfinement of nuclei and can 
be interpreted as a counterpart of the recovery of the dif- 
ferent nuclei components at higher energies (especially 
visible for the heaviest nuclei). As mentioned above, 
this decrease is slower for the largest magnetic field case. 
Above E ~ 4 x 10 18 eV, the secondary proton fraction 
increases again due to the interaction of nuclei with CMB 
photons, starting at the energy threshold ~ A x 4 x 10 18 
eV. Nuclei are completely disintegrated by CMB photons 
and turn into secondary nucleons of lower energies, re- 
sulting in an expected ~ 100% proton composition above 
E ~ 3 x 10 20 eV. Let us note nevertheless that these in- 
teractions (unlike those at lower energies) are not due to 
the cluster environment, as nuclei are barely confined by 
the cluster magnetic field at these energies; they would 
take place in a s imilar way in the intergalactic medium 
(see for instance lAllard et al.ll2008l ). 

These figures demonstrate that with a high (but 
realistic) magnetic field of B c ~ 30 /zG at the center of 
a cool core cluster, the heaviest nuclei hardly survive, 
and escape the structures only for energies around 
E ~ 10 20 eV. For lower magnetic fields, though, a 
reasonable amount of heavy nuclei can still survive 



the propagation inside clusters of galaxies. One can 
also note that no heavy nuclei survive at the highest 
energies (E > 10 20 ' 5 ) due to photo-disintegration on the 
CMB photons, meaning that the composition becomes 
100% protons in this region. The small fractions of 
light fragments (Z < 6) present at the highest energies 
after propagation inside the cluster will be rapidly 
photo-disintegrated after a few Mpc in the extragalactic 
medium. 

Figure [6] presents the cases of a non cool core clus- 
ter with central magnetic field B c = 1 /iG with the 
source located at the center (top panel), and of a cool 
core cluster with B c = 10 /iG and the source positioned 
at 100 kpc from the center (bottom panel). The upper 
panel should be compared to the lowest panel of Fig. 
as the scaling coefficient is the same in both cases, mean- 
ing that the magnetic field profiles differ here only in the 
core of the cluster (see Fig. [1]). The higher magnetic 
field and the enhanced baryonic density in the cool core 
case (see Fig. Q] and [3|) definitely play a role at lowest 
energies (E < 10 175 eV): the confinement in the core 
is more efficient and heavy nuclei are more depleted by 
hadronic interactions than in the non cool core case. The 
smaller coherence length in the center of the cool core 
cluster (see table [1]) also contributes to increasing the 
confinement time in the Kolmogorov diffusion regime, 
as T conf cx Z 1 / 3 E~ 1 / 3 B 1 / 3 \- 2 / 3 , with A the coherence 
length of the magnetic field. 

The fact that the two first panels of Fig. [5] appear more 
optimistic for the survival of heavy nuclei than for the 
non cool core case in spite of a larger magnetic field at 
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Fig. 8.— Same as Fig. [7] but with an AGN lifetime of £agn = 100 Myr. 



the center is due to our scaling of the magnetic field. 
The field being peaked at the center for the cool core 
case, an overall normalization of B to obtain a required 
value of B c at the center lowers also the intensity of the 
field at the rim of the core. We thus get a magnetic 
field that becomes quickly very faint as we move from 
the center, and in particular, it will be fainter than for 
the case of a cool core cluster normalised at B c = 1 /iG. 
These first two plots should actually be compared to cool 
core clusters normalised at B c = 0.3 /iG and 0.1 /iG 
respectively. 

It is very plausible that the acceleration sites of 
ultrahigh energy cosmic rays are not at the very center 
of the cluster of galaxies, but shifted of some hundreds 
of kiloparsecs. This would be espe cially the case of par- 
ticles are accelerated at hot spots iRachen fe Biermann 
(119931) or in lobes of radio galaxies (e.g.. lO'Sullivan et al. 
120091 and references therein), which form at a distance 
of so me hundreds of kilopa rsecs from the central black 
hole (|Bridle fc Perlevlll984f K The lower panel of Fig. [6] 
presents the resulting spectrum obtained if the injection 
of ultrahigh energy cosmic rays happens at 100 kpc of 
the center of the cluster. In a cool core cluster, the 
baryonic density falls quite steeply with the distance 
from the center: at 100 kpc, the average density is 
already two magnitudes lower than at the center. This 
explains why nuclei survive much better in this case, 
especially at low energy. The features in the heavy 
nuclei component we described qualitatively above are 
still present but with a much lower amplitude. The 
initial composition is in this case less modified by the 
propagation in the cluster environment. The position 
of the source will thus have a strong impact on the 
resulting composition of ultrahigh energy particles and 
an ample depletion of the heavy component can only be 
expected if the source is in the immediate vicinity of the 



cluster center. 

The effects of propagation in the extragalactic medium 
are not included in these results but they should not af- 
fect importantly the composition for nearby sources (at 
less than 100 — 200 Mpc). Indeed, the interactions with 
the CMB has already depleted the nuclei inside the clus- 
ter and hence further interactions will be negligible at 
ultrahigh energies. Some influence of the far infrared 
background might play a slight role for particles of energy 
E > 3 x 10 19 eV, diminishing the flux of heavy nuclei, 
but this should also be a minor effect compared to inter- 
actions experienced in the cluster. Hence at the highest 
energies, the composition should not change, we however 
expect a change in the shape of the spectrum due to the 
GZK cut-off for the remaining protons. Magnetic horizon 
effects due to the diffusion of low energy cosmic rays in 
the extragalactic magnetic field should introduce a cut- 
off for energies E < 10 17 5 eV if the source is located at 
more than ^10 Mpc, for some magnetic field configura- 
tion and strengths (seelLemo ine 2005: IKotera fc Lemoind 
I2008at iGlobus et aLll2008fl . The extragalactic magnetic 
field structure and intensity being very poorly known, we 
prefer not to introduce a large amount of new parame- 
ters, and do not compute this effect. 

In order to compare the produced cosmic ray fluxes 
with the available observed datasets of AGASA, HiRcs 
and the Pierre Auger Observatory, we computed the 
diffuse flux assuming the same mass for all clusters, 
and a number density of n s = 10 -5 Mpc -3 . Note that 
massive clusters with M > 10 15 Mq are rare objects of 
number density 2 x 10~ 6 Mpc -3 , but the density can be 
a few times 10 -5 Mpc -3 when one incl udes relatively 
small er clusters of M > 5 X 10 14 M (|Jenkins et al.1 
|200ll). Assuming that some fraction only of those 
clusters host FRI type objects, the density that is of 
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interest to us falls around ~ 1CP 5 Mpc~ 3 . We propagate 
the cosmic rays exiting the cluster environment through 
the extragalactic medium (including all the relevant 
energy loss processes), and find that the diffuse cosmic 
ray flux fits very well the observed fluxes for the range 
of parameters we chose (namely L CI = 10 45 erg s" 1 , 
which is an upper limit for the luminosity of Seyfert 
and FRI galaxies, and a 2.3 spectral index). This 
implies that higher luminosity sources or higher galaxy 
cluster densities would result in an overproduction of 
the ultrahigh energy cosmic ray flux. The normalization 
of the spectrum indeed scales as L cr x n s . One may also 
note that a harder injection spectral index will have the 
same effect of overproducing particles as compared to 
data. 

We have considered up to now spectra calculated in the 
case of a permanent regime, assuming an infinite AGN 
lifetime. It is yet of common knowledge that AGN re- 
main active only during a limited time. Studies indi- 
cate that the typical lifetime of radio sources is of order 
~ 10 7-8 y rs for FRI type galaxies and somewh at shorter 
for FRII (jParma et al.ll200l IBird et al.ll200l) . Such a 
short injection time can have a considerable effect on 
particle spectra, as one can notice that these times are 
much shorter than the confinement times of some low 
energy nuclei. 

Figures [7] and [5] present the evolution of the cosmic 
ray spectra in time, assuming a limited AGN lifetime of 
^agn = 10 Myr and 100 Myr respectively. The cosmic 
ray afterglow observed after the extinction of the source 
is due to the confinement times of different species at dif- 
ferent energies and to their variance around their mean 
value. While time goes, we observe the progressive ap- 
parition of low energy particles and of heavier nuclei. 
The variance a con f around the confinement time t con { is 
globally proportional to the latter, meaning that high en- 
ergy light nuclei have a small variance. For this reason, 
high energy protons and Helium quickly disappear as the 
source dies, while heavy nuclei with larger confinement 
time and variance remain present a much longer time. 
This leads interestingly to a heavy composition at the 
highest energies for times greater than ~ £agn + ^csc, 
where t csc is the escaping time of protons propagating 
rectilinearly from the cluster. After ~ IO^agn, the flux 
is considerably diminished at all energies. 

Such effects could be detectable if a few nearby clus- 
ters of galaxies with an extinguished AGN contribute 
significantly to the observed diffuse flux of ultrahigh en- 
ergy cosmic rays. It is also interesting to notice that 
such effects should be present in the case of any other 
type of non stationary sources embedded in magnetised 
media and injecting a mixed composition. One might 
also relate these time dependent fluxes to the absence of 
powerful sources in the arrival directions of the observed 
highest energy events: magnetised clusters of galaxies 
hosting an extinguished AGN can be emitting cosmic ray 
afterglows, and contribute to the overall observed spec- 
trum, provided that their density and source luminosity 
are high enough. 

These time dependent effects will occur according to 
the duration of the duty cycle of AGN. This value will 
set the average cosmic ray emission state at which one 
expects to observe a cluster of galaxy hosting an AGN. 



We may also notice that in the presence of more than 
one AGN at the center of the cluster, there might be 
a spread in the effective injection duration, which can 
mimic a permanent regime. In such a case, the obtained 
fluxes are those calculated in Fig. [5] and 

4.2. Secondary neutrinos 

Another signature of the propagation of ultrahigh en- 
ergy nuclei in clusters of galaxies can be found in the 
produced secondary neutrinos. One should note how- 
ever, following the previous discussion, that if the ratio 
between the durations of inactivity and activity phases of 
some AGN in clusters of galaxies is important, we might 
not observe neutrinos and photons simultaneously with 
ultrahigh energy cosmic rays. 

Figure O presents our calculated neutrino flux from a 
single cool core cluster of galaxies located at 100 Mpc, 
and embedding a central source of cosmic ray luminos- 
ity L CT = 10 45 erg s _1 that injects a mixed compo- 
sition with spectral index 2.3 with maximum energy 
E maK — 10 20 - 5 eV. We do not inject particles with en- 
ergy less than Bmin.simu = 10 16 eV and stop propagating 
nuclei when they reach this threshold: consistently, we 
do not present the flux of neutrinos below 1 PeV. These 
low energy neutrinos would in any case be rapidly over- 
whelmed by atmospheric neutrinos and make any detec- 
tion (especially for the diffuse flux) more difficult. This 
plot enables one to measure the contribution of the vari- 
ous backgrounds on the final neutrino flux. As expected, 
the hadronic interactions play a major role at all energies, 
and the cluster infrared background contributes around 
E ~ io 17 - 5 ~ 19 eV, leading to the production of neutrinos 
of energy 10 7 ~ 8 GeV. 

Figure [TU] shows the contribution of protons, helium 
and iron to the production of neutrinos. It appears that 
the contribution of protons dominates by far those of 
other nuclei, but the different components exhibit the 
same shape. The gap between the different species is 
mainly due to the relative abundances we assumed (sim- 
ilar to the Galactic cosmic ray composition) and the dif- 
ferent contributions would be closer if all the species 
had similar abundances. Even in this case however, 
the contribution of the different species would decrease 
with mass. This is due to the fact that the energy 
of the neutrinos depend on the nucleus Lorentz factor 
E v oc r^4 oc Ea/A, which means that a neutrino pro- 
duced at energy E v requires a higher cosmic ray energy 
Ea if A is large. For the typical spectral index we as- 
sume in this study (> 2) the contribution of light nuclei 
is larger. The hierarchy between the different primary 
will of course depend on the spectral index of the source, 
hard spectra giving a smaller gap between the different 
species. The lower interaction lengths of heavy nuclei 
tend to partially counterbalance this effect as heavy nu- 
clei usually suffer more interactions potentially produc- 
ing neutrinos. 

Nuclei can produce neutrinos either by interaction of 
the parent nucleus (through hadronic interactions or pho- 
topion production) or by the interactions of the sec- 
ondary nucleons (hadronic, photopion or neutron decay 
but the contribution of the latter is very low in the energy 
range we consider) . The relative contribution of the sec- 
ondary nucleons ultimately depends on the mass of the 
primary cosmic ray, as the mean free paths for photo- 
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hadronic and hadronic interactions read Ia~ oc A 1 and 
Ia p oc A~ 2 / 3 . As a result, in the case of heavy nuclei, 
the secondary nucleons are produced very rapidly, when 
the nucleus is still in the central regions of the cluster. 
This is not the case for lighter nuclei such as Helium. 
The secondary nucleons produced in the central regions 
of the cluster will encounter denser photon and baryon 
backgrounds and are therefore more likely to produce 
neutrinos. The contribution of the secondary nucleons 
to the neutrino fluxes produced inside the cluster is then 
higher for the heaviest nuclei. 

The competition between hadronic and photonic 
interaction is also an interesting aspect of the neutrino 
production in galaxy clusters. As mentioned before 
hadronic interactions dominate over photo-interaction 
for the production of neutrinos (for neutrinos produced 
below E v = 10 18 eV). In the case of protons, the 
hadronic interaction rate is far higher in the cluster 
center and photo-interaction partially compensates their 
deficit due to the flatter evolution of the photon density 
as a function of the distance to the cluster center. Let 
us note that below E — 10 16 eV, photo-interaction 
ceases to contribute to the neutrino production due 
to the energy threshold of the photopion process with 
the ambient photon background. In the case of nuclei, 
the gap between the contribution of hadronic and 
photo-interaction is larger. This can appear odd at 
first thought as the rate of photo-interactions oc A 
grows faster with the mass than the hadronic inter- 
action rate oc A 2 / 3 . Yet a significant fraction of the 
photo-interactions happens through the giant dipole 
resonance process which only produces neutrinos via 
the interactions of secondary nucleons. The photopion 
production of nuclei that can produce directly neutrinos 
whenever the produced pion is not absorbed in the 
nucleus is then screened by this less efficient process. 
On the other hand, hadronic interactions of the parent 
nucleus produce neutrinos in all the cases which explains 
the lower relative contribution of photo-interactions 
to the neutrino production in the case of nuclei. One 
may also note, that due to the different evolutions with 
the distance to the cluster center of the photon and 
baryon backgrounds the relative contribution of the 
photo-interactions will be higher when the source is 
shifted from the cluster center. 

As an illustration, the sensibility of the next gen- 
eration detector KM3NeT will be of order ~ 2 x 
10~ 9 GeV s _1 cm~ 2 for a year for a point source around 
PeV energies. We may expect from these plots that the 
flux from this galaxy cluster is typically below the cur- 
rent observable limits of neutrinos for a spectral index of 
2.3. The flux is boosted of one order of magnitude for an 
index of 2.1 and is closer to experimental sensitivities (see 
Fig. [H]), assuming the same integrated luminosity above 
1 GeV. Yet, as pointed out in the previous subsection, 
such a hard injection spectrum would overproduce ul- 
trahigh energy cosmic rays as compared to the observed 
data, by a factor of ~ 50. This means that this kind of 
central sources with high luminosities and hard spectral 
indices, more favorable for the detection of point sources, 
should have a low density if one does not want to vio- 
late the limits imposed by the total cosmic ray flux. One 



can infer from these results that detections from a single 
cluster seem difficult for the moment, but there could be 
some chance of observing diffuse fluxes. 

In this token, we calculate the expected cumulative 
neutrino fluxes for various parameter sets, for a galaxy 
cluster density n s = 1CP 5 Mpc~ 3 , assuming the same 
mass for all galaxy clusters, as we did previously for 
the cosmic ray diffuse flux. Our results are presented 
in Fig. [TTJ We do not take into account the evolu- 
tion of the source density and luminosity with redshift 
which should have a negligibl e effect for redshifts z < 0.2 
(jColafrancesco fc Blasilfl998| ). For higher redshifts, one 
should also consider the evolution of the cluster itself 
(density, infrared background, and most of all magnetic 
field), which is an impossible task seen our poor knowl- 
edge on the origin and evolution of the extragalactic mag- 
netic field. Our calculations enables us to capture the 
essential features due to key parameters, and to notice 
that all our fluxes lie around the observable threshold 
of current and upcoming experiments. Indeed, the dif- 
ferential sensibility of IceCube for diffuse fluxes is of or- 
der 1.5 x 10 -8 GeV s _1 cm~ 2 sr _1 for one year, which 
leaves room for a positive detection of signals coming 
from clusters of galaxies around 1 PeV. In this energy 
range, our fluxes are above the expected cosmogenic neu- 
trino fluxes because of the magnetic confinement and en- 
hanced baryon and photon backgrounds in the cluster 
environment. 

On the contrary for ultrahigh energies, cosmic rays are 
not significantly confined, hence we find that the neu- 
trinos produced inside the galaxy cluster by interactions 
with CMB photons (for which we take into account the 
cosmological evolution) only represent a fraction of the 
total cosmogenic neutrinos and will thus be dominated 
by them. 

It is interesting to notice the importance of the mag- 
netic confinement for the production of secondary neu- 
trinos, as well as the differences between cool core and 
non cool core clusters. As we pointed out in the previous 
section, one should only compare the non cool core case 
presented here with the cool core case at B c — 30 /iG. It 
appears that the the presence of the magnetic field en- 
hances the neutrino production of an order of magnitude, 
but there is only a slight difference between the various 
magnetic intensity and configurations. We note that a 
pure proton composition leads to a similar neutrino flux 
as compared to a mixed Galactic composition. The flux 
is actually even slightly higher, as protons produce more 
neutrinos than nuclei (see above) for the choices of as- 
trophysical parameters we have made. In the case of 
a pure proton composition, steeper source spectral in- 
dexes would be required to fit the cosmic ray spectrum, 
which means that extremely large luminosities for single 
source could be allowed (even too large to be realistic if 
one assumes a single power law down to 10 9 eV). This 
argument on the luminosities is usually alleviated by in- 
voking a change in the injected spec tru m at some en- 
ergy ( see [BcTCzinsk v fc Gazi zov 20071 and iMurase et al.l 
2008a in the context of galaxy clusters). As the fluxes 
again scale with L cr x n s , an increase in one of these pa- 
rameters could enhance the neutrino rate. One should 
yet remember that the cosmic ray fluxes would then be 
overproduced as compared to the observed data, as we 
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Fig. 9. — Neutrino flux for a cool core cluster with central mag- 
netic field B c = 10 /iG and a source located at the center, injecting 



a mixed galactic composition with L c 
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s \ at a dis- 



tance of 100 Mpc, for spectral index 2.3 (thick black) and 2.1 (thin 
violet). The contribution of the different backgrounds are indi- 
cated for spectral index 2.3: interactions with the CMB photons 
(blue dotted), with the infrared photons produced by galaxies of 
the cluster (red dashed), with the diffuse extragalactic background 
(green dash dotted) , and with the baryonic background (pink dash 
dot dotted). 




Fig. 10. — Neutrino flux for a cool core cluster with central mag- 
netic field B c = 10 /iG and a source located at the center, injecting 
a mixed galactic composition with L CT = 10 45 erg s — 1 and a spec- 
tral index of 2.3, at a distance of 100 Mpc. The contribution of 
primary protons (green dashed), helium (blue dash dotted) and 
iron nuclei (red solid) are indicated. 



calculated in section 14.11 The constraints imposed by 
the total cosmic ray flux on the diffuse neutrino flux are 
indeed quite stringent. However the dilution of the flux 
due to the limited AGN lifetime can be used to justify an 
increase of the luminosity by a factor £agn /t cycle, where 
icycie is the periodic duration of an AGN cycle. 

Our r esults are consist e nt wi th the analytical treat- 
ment of iBerezinskv et all (|1997l ) - and with the study 
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Fig. 11. — Diffuse neutrino fluxes obtained with galaxy cluster 
density n B = 10~ 5 Mpc -3 and AGN cosmic ray luminosity L CT — 
10 45 erg s — 1 . A mixed composition is injected at the center of 
the non cool core cluster with B c = 1 fiG (red thick solid), and 
in cool core clusters with B c = 30 fiG (black thick solid), B c = 
10 fj,G (black thin solid), B c = 3 fiG (black dotted) and without 
magnetic field (green dash dotted). We also present the cases of 
a pure proton injection at the center (blue long dashed) and a 
mixed composition injected at 100 kpc from the center of a cool 
core cluster of B c = 10 /iG (pink solid). 



of iMurase et al.l (|2008aD though they assumed different 
physical parameters. A rough order of estimate on the 
neutrino flux J u around PeV energies in the case of a 
pure proton composition, assuming that the hadronic in- 
teractio ns are the dominant interaction process can be 
written (|Murase et al.ll2008af) : 



E 2 J V {E) 
fpp 



-2 



< 0.7 x KT 11 GeV s" 1 cm" 

rcr \ 
E16 ^ , (1) 



D 



-2 



-1 



E.16 



\2A x 1CT 3 / V 100 M P<V \10 43 ergs 

where D is the distance to the source, L 
10 43 erg s" 1 the cosmic ray luminosity at E = 
10 16 eV (corresponding roughly to a value of 16 = 
10 45 erg s _1 for a minimum injection energy of E m i n — 
10 9 eV, with spectral index 2.3) and / pp the ef- 
fective optical depth for the proton-proton interac- 
tions at energy E ~ 10 16 eV. This latter quan- 



tity can be written: f p 



0.8 (TppTlNCte 



2.4 x 



10~ 3 (n#/10~ 4 - 5 cm~ 3 )(£ CS c/l Gyr), assuming a constant 
baryonic density nn and escape time t csc throughout the 
cluster. 

Our fluxes are l ower than those calculated by 
Ide Marco et~aH (|2OO60 in the energy range range between 
10 16 and 10 18 eV, and show an overall difference in the 
shape of the energy spectrum. This discrepancy stems 
mainly, as already mentioned in section 12.31 from their 
choice of very bright infrared galaxy SED (instead of 
elliptical galaxy in the present study) to calculate the 
cluster photon background. Furthermor e, hadronic in- 
teract ions were not taken into account bv lde Marco et aD 
(|200l) . 

The neutrino fluxes presented in Fig.s [51 [TO] and QT] 
do not take into account the limited AGN lifetime and 
assume a permanent emission regime. This is justified 
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for the highest energy cosmic rays that produce neutri- 
nos through interactions with the CMB photons and that 
are not trapped inside the cluster: neutrino production in 
this case should thus happen quickly after the injection. 
We checked that it is also the case for the relatively lower 
energy particles. Indeed, most of the PeV energy neu- 
trino flux is produced in the central region of the cluster 
shortly after injection. 

4.3. Secondary gamma rays 

Secondary gamma rays can also be a signature of the 
propagation of ultrahigh energy protons or nuclei in clus- 
ters of galaxies. As for the neutrinos, the simultaneous 
observation of charged particles and of gamma ray pho- 
tons from a cluster will depend on the duration of the 
life cycle of the source. 

Very high energy charged and neutral pions are pro- 
duced via hadronic and photo-hadronic interaction pro- 
cesses. Neutral pions decay into high energy gamma 
rays, while charged pions lead to high energy electrons 
and positrons that lose their energy through synchrotron 
and Inverse Compton emissions. Generated gamma rays 
are attenuated by the pair-creation process on the inter- 
galactic photonic background during their propagation, 
while electrons and positrons typically lose most of their 
energies via synchrotron and Inverse Compton losses lo- 
cally and sufficiently low energy pairs are trapped by the 
magnetic fields of clusters. 

In our propagation code, the gamma ray, electron and 
positron production via pion decay is t reated in a dis- 
crete manner, as in lAllard et all (|2006fl for nuclei pro- 
jectiles with A > 1 and with SOPHIA for proton and 
neutron cosmic rays. Electron and positron pair cre- 
ation throu£h_rjhoto^adronic processes is implemented 
as in lArmeng aud et al.l (|2006f ). For each time step (cho- 
sen to be much larger than the typical mean free path of 
pair photo- production processes) , we assume that an en- 
semble of electron and positron pairs are generated with 
energies distributed according to a power law, and max- 
imum energy depending on the primary particle. As we 
are interested in the cascad ed photons, the use of an im - 
proved pair spectrum as in Ke lner fc Aharonianl (2008) 
is not expected to change our results. 

Electromagnetic cascades are treated as post-analysis, 
using the radial distribution of electrons, positrons and 
photons produced in the simulation, inside the cluster 
of galaxies. Very high energy photons propagate in the 
direction of ultrahigh energy cosmic rays, while lower en- 
ergy photons are radiated almost isotropically. Here, for 
simplicity, we approximate that photons propagate in the 
radial direction. Although this is a crude approximation, 
it is justified in so far as the energy attenuation length 
due to the CMB (at ~ PeV) is much shorter than the 
cluster size of ~ Mpc, and it is enough to obtain qualita- 
tively appropriate results (Asano 2009, private communi- 
cation). The resulting high energy gamma ray fluxes are 
calculated by so l ving kinetic e quations (see for example 
Aharonian 2002; Murase 2009 and references there in). 

In Figure I12[ we show the spectra of produced gamma 
rays from a cool core cluster with B c = 3 /iG, which 
is located at a distance of 100 Mpc, with a source of 
luminosity L CI = 10 45 erg s _1 and of spectral index of 
2.3 (the parameters here are identical to those used for 
neutrino fluxes in Fig. [9]). The case of a non cool core 



cluster with B c = 1 /iG is also shown in Fig. [T3] Only 
the contribution of cosmic rays with E > 10 16 eV is 
shown in these plots. Very high energy electrons mainly 
radiate their energy through synchrotron emission and 
lead to the production of gamma rays around keV-GeV 
energies. On the other hand, gamma rays produced via 
Inverse Compton emission are expected above the TeV 
range, which is indeed observed in the figures (thin red 
line in Fig. [12] and green dashed lines in Fig. [T3|). 

Thick red line in Figure [T^] and red lines in Fig- 
ure [T3] indicate the flux obtained after propagation of 
photons, electrons and positrons in the extra-cluster 
medium. Gamma rays with > PeV gamma rays are sig- 
nificantly attenuated by the cluster photon field and the 
CMB (where the cluster photon field itself causes only 
a minor effect). Above ~ 10 TeV, photons also suffer 
from interactions with the diffuse infrared, microwave 
and radio backgrounds. Thus, they are largely absorbed 
and cascaded into lower energy gamma rays especially 
for high redshift clusters of galaxies. In our calcula- 
tions, we set the intergalactic magnetic field strength to 
Big = 10 -13 G, which allows us to consider the cascades 
as unidimensional as long as the coherent length is short 
enough so that the deflection angle is small, 9b -C 1 
(|Murase et al.l l2008cf) . For stronger magnetic fields or 
longer coherent lengths, one must consider the deflec- 
tions experienced by electrons and positrons, which will 
lead to a dilution of the flux. The synchrotron emis- 
sion from electrons and positrons near the source (green 
dashed line in Fig. [T2"]) as well as photons that are di- 
rectly produced by cosmic rays below some TeV will not 
be affected, as the Universe appears transparent to those 
particles. It should be noted however that the contribu- 
tion of photons produced by cascades in the intergalactic 
medium is relatively limited below ~ 1 TeV as compared 
to the flux produced in the cluster itself (see Fig. [T2"]) . 

We also assume that our clusters of galaxies harbour a 
steady AGN as a central source. Note that considering 
a finite AGN activity time could lead to a diminished 
gamma ray flux around TeV energies, if the intergalactic 
magnetic field is strong enough to spread the cascaded 
gamma ray emission over a much longer time than the 
AGN lifetime. These issues will eventually depend on the 
rate of magnetic enrichment of the intergalactic medium. 
As for the synchrotron emission directly produced in the 
cluster, and that constitutes the major gamma ray con- 
tribution at lower energies, a finite AGN lifetime could 
lead to a complete temporal decorrelation of the gamma 
ray flux from the cosmic ray flux, according to the ratio 
between the source activity and inactivity durations. 

We observe in Fig. [12] and [13] the same effects as for 
neutrinos, depending on the configuration and the inten- 
sity of the magnetic fields. For gamma rays, one may no- 
tice an additional contribution of the intensity of the field 
on the synchrotron emission: for stronger fields, elec- 
trons and positrons produce a larger amount of photons 
of higher energy (the energy loss length of an electron of 
energy E e is indeed x syn oc B~ 2 E~ X and the synchrotron 
emission peaks at energy E^ syn oc BE%). These features 
are not striking in figure [T3l where the magnetic field is 
always strong enough to have an important synchrotron 
emission. The scenario in which the cluster magnetic 
field is set to zero suppresses completely the flux below 
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Fig. 12. — Gamma ray emission from a source placed in a cool 
core galaxy cluster with B c = 3 fiG, a luminosity of L cr = 
10 45 erg s — 1 , with spectral index 2.3, and located at a distance of 
100 Mpc. Blue dotted lines represent electron and positron (thin) 
and photon (thick) fluxes produced in the cluster by cosmic rays of 
energy E > E m i n 8 j mu = 10 16 eV. Green dashed lines indicate the 
contribution of synchrotron emission that occurred in the cluster. 
The red thin line is the gamma ray flux obtained after propagation 
in the cluster of galaxies, that one would observe at a distance of 
100 Mpc in absence of electromagnetic cascades in the extra-cluster 
medium. The thick red line indicates the photon flux observed at 
100 Mpc after electromagnetic cascades, assuming an extra-cluster 
average magnetic field intensity of Big = 10 — 13 G. 
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Fig. 13. — Comparison between the fluxes obtained for a cool core 
cluster with B c = 3 fiG (thick lines) and a non cool core cluster 
with B c = 1 fi G (t hin lines) - the other parameters are identical 
to those of Fig. 1121 Green lines represent the flux obtained after 
propagation in the cluster, as one would observe it at 100 Mpc in 
absence of electromagnetic cascades in the extra-cluster medium. 
Red lines give the gamma ray flux observed at 100 Mpc after elec- 
tromagnetic cascades, assuming an extra-cluster average magnetic 
field intensity of Bi G = 10" 13 G. 



~ 10 GeV. 

We examine the contribution of cosmic rays of en- 
ergy E < i? m in,simu = 10 16 eV to the flux of gamma 
rays in the GeV-TeV energy range. The propagation of 
these low energy particles has not been calculated nu- 
merically because of the excessive calculation time that 
it implies. It is however possible to calculate the spa- 
tial distribution and the energy spectrum of these cos- 
mic rays, assuming that they are in a diffusive regime, 
which is appropriate in this energy range. As gamma 
rays are mainly produced by protons through hadronic 
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Fig. 14. — Contribution to the gamma ray flux of low energy 



cosmic rays with E < E n 



10 16 eV, for a cool core cluster 
with B c = 3 /iG — the other parameters are identical to those of 
Fig. 1121 Thick lines represent the calculated emission for cosmic 
rays of energy E > E m i n ,simu, and the thin lines the emission in- 
cluding low energy cosmic rays modeled semi-analytically in the 
diffusive approximation. The cosmic ray injection time is set to 
3 Gyr, assuming that the sources emit particles almost continu- 
ously. 
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FiG. 15. — For illustrative purposes: gamma ray fluxes produced 
by low energy cosmic rays (E < -E m i n ,simu) assuming spectral 
indices of 2.1 (red solid) and 2.3, with different minimum injec- 
tion energy _E min = 1, 10, 100, 10 3 GeV. The source luminosity is 
assumed to be L cr = 10 45 erg s — 1 for spectral index 2.3, and 
L cr = 2 X 10 43 erg s — 1 for spectral index 2.1 between 1 GeV and 
10 20 - 5 eV. 



interactions in this energy range (the discussion is sim- 
ilar to the case of neutrinos, see previous section), we 
consider only proton-proton interactions in our calcula- 
tions. We assume an initial Gaussian distribution that 
extends over ~ 1 kpc in the cluster, and we further sup- 
pose that the central source injects particles at a rate that 
maintains this gaussian distribution at all times. The 
diffusion equation is solved by using a Crank-Nicholson 
scheme, adopting the following diffu sion coeffici e nt com - 
puted according to the results of iCasse et al.1 |2002) : 
D = 1.2 r L c [r L // c + 0.1(r L /^ c ) ~ 2 / 3 ] , where r L is a Lar- 
mor radius and l c is a coherenc e length of the magnetic 
field (|Kotera fc Lemoind[2008H Appendix A). 

In the inner region of the cluster (~ 10 kpc), we ob- 
tain an energy spectrum steeper than the injected one, as 
low energy particles are accumulated while higher energy 
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particles escape from the inner region earlier. The op- 
posite effect is observed in the outer regions (~ 1 Mpc), 
where low energy particles are less numerous. One can 
note that the calculated initial photon spectrum slope is 
of order —2.6, which corresponds as expected to the com- 
bination of the injection slope (—2.3) and to the evolution 
of the confinement time with energy (—1/3). We evalu- 
ate the gamma ray flux produced by proton-proton inter- 
action s thro ugh the formulations o f iPfrommer fc"li lnfilin 
(|2004[ ) and lAndo fe Nagal ([20081 ). Finally, we match 
the resulting power-law spectrum with the numerically 
calculated one by interpolating the value between E = 

8 x 10 15 eV (corresponding to the maximum value of 
our calculation under the diffusion approximation) and 
E = 10 16 eV (corresponding to the minimum value of 
the higher energy spectrum). Such a rough matching is 
sufficient for demonstrative purpose. 

Our results are shown in Fig. [14] and 1151 which show 
that gamma ray signatures of ultrahigh energy cosmic 
rays are overwhelmed by gamma rays produced by 
lower energy cosmic rays. The situation looks less 
critical for harder cosmic ray spectra, but the photons 
emitted by ultrahigh energy cosmic rays will still 
be one or two orders of magnitudes below the low en- 
ergy particle contribution for spectral indices of 2.0 — 2.1. 

We discuss in what follows the detectability of pho- 
tons from clusters of galaxies. The future Cherenkov 
Telescope Array (CTA) is expected to have a sensitivity 
of order ~ 10~ n GeV s _1 cm~ 2 for 100 hours for a point 
source around TeV energies. For an extended source, the 
sensitivity becomes weaker by a factor of 0/Opsf, where 

9 and $psf are the angular extensions of the source and 
of the point spread function respectively. The latter be- 
ing of order of a fraction of arcminute around 1 TeV, 
we might assume that a typical cluster of galaxy ob- 
served at a distance of 100 Mpc leads to an angular size 
of 6 ~ 10 6>psf — 100 6*psf, which corresponds roughly 
to a sensitivity of - lO" 10 - 10~ 9 GeV s^cm" 2 . The 
Fermi satellite has a lower angular resolution (of a frac- 
tion of degree) around GeV energies, leading to a sensi- 
tivity reaching ~2 x 10~ 10 GeV s _1 cm~ 2 for a year for 
a typical cluster of galaxies at 100 Mpc in this energy 
range. 

Figure [T^] shows that the gamma ray signal coming 
from ultrahigh energy cosmic rays would be hardly de- 
tectable by any of these instruments. This conclusion 
also applies to the future detectors HESS-2 and MAGIC- 
2. Those signals might moreover be overwhelmed by 
a stronger contribution of low energy cosmic rays, as 
pointed out in Fig. [T3J 

The gamma ray flux from cosmic rays of low energy 
might be detected by the experiments cited above, but 
they should first be compared to existing observational 
limits. The EGRET point source limit is of order 3 — 
4 x 10~ 9 GeV s" 1 cm" 2 at 100 MeV, for a cluster lik e 
Coma located at ~ 100 Mpc (jPfrommer fe Enfilin| [2004) . 
We thus note that the values presented in Fig. Qj)] for 
spectral index 2.3 violate this limit in the energy range 
~ 10 8-10 eV. Besides, an estimate of the diffuse gamma 
ray flux using the same method as for neutrinos leads to 
a value of ~ 5 x 10~ 5 GeV s _1 cm~ 2 sr _1 at 1 GeV for 
a spectral index of 2.3, which is also above the diffuse 



EGRET limit ([Osborne et al.|[l99"l . 

It might be reminded however, that these plots were 
made under the assumption of a stationary cosmic in- 
jection regime, and that considering a limited AGN life- 
time and particular configurations of sources may lower 
enough the flux to avoid the violation of the EGRET 
limit. These problems can also be overcome in the case 
of a harder cosmic ray injection spectrum (< 2.1), that 
can be either a continuous power law extending to ul- 
trahigh energies (in which case one should remember 
that the luminosity or the density of the source should 
be lowered by a factor ~ 50 in order not to overshoot 
the observed cosmic ray data, implying a luminosity of 
L CT = 2 x 10 43 erg s" 1 between E = 10 9 - 10 20 5 eV), or 
a broken power law. In the case of a 2.1 low energy injec- 
tion spectrum, our calculated diffuse flux from clusters of 
galaxies would contribute to ~ 20% of the EGRET dif- 
fuse flux. Such a flux should be also detectable by Fermi 
around GeV energies. In any case, future observations 
will be able to put constraints on our predictions. 

We also made for these calculations the assumption 
that AGN in clusters of galaxies were the major con- 
tributors to the total ultrahigh energy cosmic ray flux 
in order to set the source luminosity and density. One 
may however alleviate this hypothesis and invoke a lower 
source density or luminosity leading to a better consis- 
tency with the observed limits. Moreover, these results 
raise the question of the minimum energy of cosmic rays 
injected in the cluster of galaxies. It is indeed not obvi- 
ous that cosmic rays can escape the acceleration site of 
AGN at an energy as low as 1 GeV, due for example to 
high magnetic con finement times and sh orter adiabatic 
energy loss times ([Berczins kv et al1 l2006). If the cosmic 
ray injection starts at an energy greater than 100 GeV 
(as presented in Fig. [T5"|) , the fluxes could be compatible 
with observed data, even with spectral indices of 2.3, but 
their detectability by current and upcoming instruments 
might be compromised. 

Ultrahigh energy photons could al so be a pote ntial 
probe of ultrahigh energy cosmic rays (Murasc 200£j) and 
they could even produce a pair halo after they leave a 
cluster. However, our results imply that their detection 
by current and future detectors such as CTA does not 
seem easy as well, unless the source is rather close or 
luminous. 

In this paper, we did not evaluate the energy spectra 
of very low energy photons, especially in the radio 
band. Values for the diffuse radio emission flux from 
some cl usters of galaxies were r e ported by several au- 
thors (iGiovannini et all 119931 : iGiovannini fc Ferettil 
1 20001 IFeretti et al l 120041 : iBagchi et al I 120061 : 
iFeretti fc Giovanninil 120081 : IBrunetti et all [2008). 
In a few years, the Low Frequency Array (LOFAR) and 
the Long Wave length Array (LWA) will even observe 
clusters at low radio frequencies, and possibly discover 
the bulk of t he cluster scale synchrotron emission in 
the U niverse ([EnBlin fc Rottgerind [20021: IBrunetti et~aTI 

Hoop). 

Finally, it should be noted that these low energy 
photons as well as < TeV gam ma rays may come 
from leptonic components (e.g.. lEnsslin fc Biermannl 
1991: lAtovan fc V51kl [2000t iTotani fc Kitavamal 12001 
Kushnir fc Waxmanll2009l . and references therein) rather 
than hadronic ones. One promising way of distinguish- 
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ing them would be to observe the spatial distribution of 
gamma rays, since hadronic gamma rays trace th e gas 
while leptonic ones trace the shocks (|Mi niati 2003|). 

5. CONCLUSION 

We studied the propagation of nuclei in a magnetised 
cluster of galaxies using a complete propagation code and 
based on a detailed study of the physical properties of 
clusters of galaxies. Due to their strong magnetic fields 
and dense matter and radiation backgrounds in the cen- 
tral regions, galaxy clusters are in principle hostile envi- 
ronments for cosmic ray nuclei. It is found that the nuclei 
survival depends mostly on the strength and the profile 
of the magnetic field, as well as on the position of the 
source. Indeed, heavy nuclei can be efficiently suppressed 
when the source is found at the very center of a cool core 
cluster for the strongest magnetic field normalizations. 
Light and intermediate nuclei are found to be much less 
affected by the confinement and the dense central pho- 
ton and baryon backgrounds even when strong magnetic 
fields are at play. In the case of lower central magnetic 
fields or if the source is shifted from the central regions, 
the survival of the heaviest elements becomes easier. 

We also examined the effect of the limited lifetime of 
the central source. The resulting spectrum from a single 
source is found to be time dependent due to the rigidity 
dependence of the confinement time in the magnetized 
regions. As a consequence, at a given energy, the com- 
position becomes heavier, the flux fainter with time and 
the overall spectrum gets softer. These effects, though 
expected for any source surrounded by a magnetized re- 
gion, should not affect the overall cosmic ray spectrum 
and composition providing that a sufficiently large num- 
ber of source is present within the cosmic ray horizon at 
a given energy. Let us note that depending on the dura- 
tion of the life cycle, one can observe preferentially the 
contribution of both cosmic rays and secondary neutri- 
nos and gamma rays, or the former without the latter, 
or even the latter without the former. 

Concerning secondary neutrinos, we have shown that 
their production can be significantly enhanced by the 
confinement of the cosmic rays (below ~ Z x 10 18 eV) 
in the central regions of the cluster. As a consequence, 
neutrinos below 10 17 eV, produced mainly via hadronic 
interactions or interactions with ambient infrared and 
optical photons, could form a detectable diffuse flux for 
IceCube or KM3Net if a significant fraction of all the 
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UHE cosmic rays originates from the central regions of 
galaxy clusters. At the highest energies however, cosmic 
rays are barely confined by the cluster magnetic fields 
and as a result, the diffuse neutrino flux from galaxy clus- 
ter above 10 18 eV should only represent a small fraction 
of the whole diffuse cosmogenic neutrino flux, preventing 
the signature of their origin. The flux from single sources 
(except if they are exceptionally bright) should be well 
below the current observation capabilities. 

Gamma rays are another potential signature of the 
presence of cosmic ray accelerators in Galaxy clusters. 
We find that the fluxes produced by very and ultrahigh 
energy cosmic rays should remain typically below the 
reach of current and planed experiments for reasonable 
assumptions on the sources luminosity, and correspond- 
ing source densities of n s > 1CF 6 Mpc -3 . Moreover, 
due to the combined effects of the cascading of the pro- 
duced photons and the evolution of the confinement time 
with the rigidity (coupled to the flatness of the hadronic 
interaction cross sections), this contribution is likely to 
be overwhelmed by the one of lower energy cosmic rays. 
If the latter component is released in the intra-cluster 
medium, ample gamma ray fluxes, testable by satellites 
or ground based experiments could be produced. 

As a conclusion, multi-messenger observations could 
be used to test the presence of cosmic ray accelerators 
in galaxy clusters. Though ultrahigh energy cosmic rays 
are more likely to be probed by their arrival directions, 
very high and low energy cosmic rays could be revealed 
using neutrino and gamma ray observations. 
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